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ABSTRACT 

We have developed a new numerical approach to the dynamics of minor 
bodies and dust particles, which enables us to increase, without using a 
supercomputer, the number of employed particle positions in each model up to 
10 10 — 10 u , a factor of 10 6 — 10 7 higher than existing numerical simulations. 

We apply this powerful approach to the high-resolution modeling of the 
structure and emission of circumstellar dust disks, incorporating all relevant 
physical processes. In this Letter , we examine the resonant structure of a 
dusty disk induced by the presence of one planet of mass in the range of 
(5 • 10~ 5 — 5 • 10 -3 ) A/*. It is shown that the planet, via resonances and 
gravitational scattering, produces (i) a central cavity void of dust; (ii) a trailing 
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(sometimes leading) off-center cavity; and (iii) an asymmetric resonant dust 
belt with one, two, or more clumps. These features can serve as indicators of 
planet (s) embedded in the circumstellar dust disk and, moreover, can be used 
to determine the mass of the planet and even some of its orbital parameters. 
The results of our study reveal a remarkable similarity with various types of 
highly asymmetric circumstellar disks observed with the JCMT around Epsilon 
Eridani and Vega. 


Subject headings: circumstellar matter-dust-planetary systems-stars: individual (e 
Eridani, Vega)- 
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1. INTRODUCTION 


1.1. Current approaches to the modeling of dusty disks and their limitations 


Dusty disks around stars are a very common cosmic phenomenon (Beckwith k Sargent, 
1993; Marcy, Cochran k Mayor, 1999; Woolf k Angel, 1998). Both the number and variety 
of circumstellar dusty disks being discovered by the HST and other telescopes grows swiftly, 
bringing new data virtually every month. These disks, with scales of tens to hundreds (up 
to a thousand) of AU, possess a great variety in structure, sometimes they have a central 
‘hole’ void of gas and djast, and often are highly asymmetric. Those disks are thought 
to accompany planetary systems and, at the same time, to hide them from the observer. 
Knowledge of dust characteristics is of prime importance for future NASA’s space missions, 
such as the NGST (Next Generation Space Telescope - see Mather et al. 1997) and the 
Terrestrial Planet Finder (TPF, 1999). It is generally expected that a significant limitation 
to unambiguous planet detection and study will be contaminating thermal emission from 
dust in the target systems. We propose to turn this hazard into an advantage for detecting 
and characterizing planetary systems. 

The development of a detailed theoretical model of exo-zodiacal dust clouds is 
complicated by 

• the wide variety of relevant physical processes, such as (i) dissipative effects (Poynting- 
Robertson drag, stellar wind, and aerodynamic drag), (ii) resonant interaction with 
planets, (iii) gravitational scattering by planets, (iv) evaporation and sputtering of 
dust particles, (v) mutual collisions in the cometary and dust populations, and 

• uncertainties in the distribution of comets or other bodies as the major sources of 
dust in planetary systems. 

Because of the complexity of the problem, the numerical approach is usually taken in 
modeling the IPD. Up until now, numerical models suffered from the limited number of 
particles that could be used in the computations. For instance, the first results in modeling 
dynamics of dust from the Kuiper belt were obtained by Liou, Zook, k Dermott (1996) with 
just 80 particles. Even now. numerical models treat several thousands of particles at most. 
While helpful in studying the dynamics of individual dust particles, these models cannot 
easily compute the large-scale structure of an IPD cloud , maps of the zodiacal emission of a 
circumstellar disk, etc. 
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1.2. Present Numerical Approach 

Recently, we have proposed a novel, very effective approach to numerical modeling 
of distributions of test particles in an external gravitational field (Ozernoy, Gorkavyi, & 
Taidakova, 1998 = OGT 1998), Gorkavyi, Ozernoy, Mather, & Taidakova = GOMT 1999). 
In our approach, the particle-number limitation is solved, which is decisive to provide 
reliable numerical simulations. In brief, our approach (which has a number of common 
elements with the ‘particle-in-cell’ computational method) is as follows: Let us consider, for 
simplicity, a stationary (in the frame co-rotating with the planet) particle distribution. The 
locus of the given dust particle’s positions (taken, say, as 6 • 10 3 positions every revolution 
about the star) are recorded and considered as the positions of many other particles 
produced by the same source of dust but at a different time. After this particle ‘dies’ (as 
a result of infall or ejection from the system by a planet-perturber), its recorded positions 
sampled over its lifetime form a stationary distribution as if it were produced by many 
particles. Typically, each run includes 10 4 - 10 5 revolutions, i.e. ~ 10 8 positions of a dust 
particle, which is equivalent, for a stationary distribution, to 10 8 particles. If we allow for 
100 sources of dust (we can include even a larger number of sources), we deal, after 100 
runs, with ~ 10 particle positions as if they were real particles. This approach brings an 
factor of 10 6 improvement in the detail of a model! In the present Letter, we will not keep 
information about the dynamical path of each particle (as we did in OGT 1998 and GOMT 
1998e), but instead, we will immediately sort this information into 8 • 10 6 spatial cells (each 
cell containing 10 3 - 10 4 particles), thereby forming a 3D grid that models the dust cloud 
around the star (GOMT 1999). Our approach makes it possible to study, besides stationary 
processes, certain non-stationary processes as well, e.g. evolution toward steady-state 
distributions (OGT 1998, GOMT 1999), dust production from non-steady sources, decrease 
in particle size (due to evaporation) and number (due to collisions), etc. 

In order to see the merits of our approach, we compare it, in Table 1, with that 
of Roques, Scholl. Sicardy, & Smith 1994 (hereinafter, RSSS), one of the best modeling 
of planetary dust disks so far. As can be seen, our approach offers greatly advanced 
opportunities for modeling circumstellar dust disks. In Sect. 2. we briefly duscuss the 
theoretical background. We explore the possible variety of dust disks with one embedded 
planet (Sect. 3) and then we make a preliminary interpretation of observations of certain 
circumstellar disks (Sect. 4). Our conclusions are given in in Sect. 5. 
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2. CIRCUMSTELLAR DISKS: THE ROLE OF RESONANCES 


In the Solar system, the Kuiper belt objects (‘kuiperoids’), along with comets and 
asteroids, serve as major sources of dust. These constituents are thought to be rather 
common for many circumstellar disks, which should contain not only planets but residual 
material such as kuiperoids, comets, and dust produced by these objects. Both the P-R 
and stellar wind drags (even stronger than in the Solar system whenever the star is more 
luminous) as well as residual gas in the disk tend to induce dust inflow toward the star. As 
the dust passes by the planets in its infall, it interacts with them by accumulating in the 
outer planetary resonances. Therefore, as we demonstrate below, the resonant structure in 
the dusty circumstellar disks could serve as an efficient means of planet detection. 


2.1. Density Distribution of Dust in the Disk Without a Planet 


It will be instructive to start with a limiting case of a circumnuclear disk composed 
only of sources of dust (no planets yet). We have computed a 3D model of such a disk 
produced by 25 kuiperoidal sources of dust experiencing the P-R drag. For illustrative 
purposes, the distribution of exo-kuiperods is taken similar to the kuiperoids in the Solar 
system. At distances from the star less than ca. 30 AU, the surface density of. dust is 
found to be almost constant forming a plateau down the star, whereas the number density 
increases inward as ~ 1/r. This result turns out to be in good accordance with our 
earlier analytical work (GOMT, 1997), which addressed the dust production by asteroids in 
the Solar system. Interestingly, the latter dependence differs substantially from that in the 
outer Solar system, where the presence of four giant planets results in « const (GOMT 
1999). 


2.2. Resonance Dust Rings Induced by Exo-planets 


The major factor which makes visible the re-distribution of dust in the circumstellar 
disk due to the presense of embedded planet(s) are resonances. The existence of a resonance 
ring associated with the Earth's orbit was predicted by by Jackson & Zook (1989); such 
a ring was indeed discovered in the IRAS (Dermott et al. 1994) and G’OBE (Reach et al. 
1995) data. Below, we focus on the conditions favorable to formation of resonant rings near 
the planetary orbits; this analytical consideration is very helpful in planning the numerical 
computations for a variety of environments. 
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2.2.1 Resonant capture 

Weidenschilling and Jackson (1993) found a criterion for the resonant capture of 
particles into outer planetary resonances (2D case, the planet on a circular orbit): the 
capture into a resonance becomes possible for those resonances for which the gain in 
angular momentum is equal to the angular momentum loss due to dissipation, i.e. when the 
‘capture parameter’ exceeds unity. This parameter depends upon the resonance number, 
the planetary mass, m p<) orbital radius, and a value of (3 « Z*/(A/», that is the ratio of 
the stellar light pressure and the gravitational force applied to the dust grain of radius r (in 
pm) (Z* and M. being the star’s luminosity and mass in solar units). Generally speaking, 
the ‘capture parameter’ depends also upon the eccentricity of the planet, e, and inclination 
of particles, i. For a giWrn mass of the planet and size of dust particles, there are certain 
resonances favorable for the particle’s capture: distant resonances are weaker and unable 
to intercept particles, whereas resonances close to a planet, which are stronger, do not get 
enough particles. The resonances determining the overall shape of a resonant ring and its 
asymmetry are entirely different for the Earth, which tends to capture particles into the 
resonances 6:5, ..., 11:10, etc., than for Neptune, which mainly captures particles into the 
resonances 2:1, 3:2, 4:3, etc., or Jupiter, whose major interceptors are the resonances 2:1, 
3:1, 4:1, etc. 

2.2.2 Resonant lifetime 

Another important factor influencing the shape and structure of resonant belts is 
the dynamics of resonant particles, which can be characterized by the two closely related 
parameters: the lifetime spent in a resonance and the maximum eccentricity, e max , 
attained by a resonant particle. (The eccentricity grows due to the fact that the PR-drag 
decreases the particle’s angular momentum even if the particle is stuck in a resonance, 
i.e. keeps its a = const). Besides an increase in e, the resonant particles change, in a 
complicated way, their i - cf. Liou Sc Zook (1997). Weidenschilling Sc Jackson (1993) found 
fimax = [0.4/ ( 1 + J)] 1 ^ 2 (j being the resonance number), at which point, the pericenter of 
a particle’s orbit becomes smaller than the planet’s orbit ( q % O.S5a p /) so that a further 
balance in angular momentum becomes impossible anymore and the particle leaves the 
resonance. 

Dust particles captured into a resonance can reach e ~ e max and spend an appreciable 
amount of time in this metastable state (Weidenschilling Sc Jackson 1993). Obviously, the 
appearance of this resonant belt only depends on a rea and e max . which are the same for all 
particles in the belt, regardless their initial a Q and eo- 

For the Lindblad resonances (j + 1) : j as well as for the j : l resonances, which 
are of major interest to us, the lifetime of particles in the above metastable state rapidly 



increases with the distance from the planet, whereas the probability of the resonant capture 
decreases. Hence, it would be a reasonable approximation to examine the contribution of 
various resonances into the the observed structure of the circumstellar disk by assuming 
that this structure is induced by just one of the outer resonances. Indeed, the outermost 
resonances are unable to capture particles, while the particle lifetime in the innermost 
resonances is much shorter. This approximation is employed in Sect. 4 when our simulated 
disks are compared with the observed circumstellar disks. 

2.2.3 Mass of a Resonant Dusty Structure 

The total mass of dust contained in the particular resonance is the product of the 
probability for the dust particle be captured into the resonance, the lifetime spent in 
the resonance, and the total dust flow per unit time in the vicinity of the resonance, i.e. 

Mj = WjTj(dM/dt). Having introduced, for convenience, the mass of ‘background’ dust, 
M&, which is common for all the resonances, one gets M 3 jMb ~ Wj / ( S 3 — 1) ( Tj/r PR ), where 
Sj is the probability for the dust particle be ejected from the system due to gravitational 
scattering on the planet, and rpp is the characteristic time of the particle’s semimajor axis 
change due to the Poynting-Robertson drag. 


3. MODELLING CIRCUMSTELLAR DUST DISKS 

We consider the simplest planetary system consisting of just one planet orbiting a 
star of mass M ~ A/@ and sitting at a eccentric orbit of a semi major rx i s q 30 AU. 
To explore the role of a particular resonance, we put the source of dust just outside of 
this resonance. Since the source of dust is assumed to belong to the Kuiper belt objects, 
its initial eccentricity and inclination were taken typical for kuiperoids (disk population): 
eo « 0.07 and t’o % 2 — 3°. We have computed 40 model disks in order to explore the effects 
of various resonances adopting as parameters: (i) the mass of the planet; (ii) its orbital 
eccentricity; and (iii) particle size coupled with stellar luminosity and mass, as described by 
the parameter /?, the ratio of the stellar light pressure and the gravitational force applied to 
the dust grain; 0 ss 0.3 Lj(Mr) depends upon the particle’s radius r (in jj.m) and the star’s 
luminosity L and mass M (in solar units) assuming dust grains of density 2 gem -3 . Each 
model is based on just one run that reveals the distribution of dust particles between their 
capture into a resonance and eventual leaving it. (In fact, we explored 10 to 15 runs with 
various inital conditions and found the resulting distribution to be only slightly varying 
from one run to another. This confirms the conclusion reached at Sect. 2.2.2 that the 
appearance of a resonant belt is not sensitive to «o and e<).) In those cases when the particle 
lifetime in the resonance was too long (exceeded 2 • I0 (i planet’s revolutions, the run was 
stopped. 1 his set of models was easily constructed, as the the ('PI 1 time to compute any 



one model on a PC typically did not exceed a few hours (or dozens of hours for very small 
0's as well as very distant resonances). Below, we describe some of these models, in order 
to reveal the qualitative dependencies upon the above parameters. 


3.1. Dependence of Structure Upon Particular Resonances 

How one planet embedded into a circumstellar disk modifies its dust distribution 
through resonances, is shown in Figs. 1 and 2 as a set of 12 models that can be used to see 
the effects of different parameters. In each panel, the surface density of the disk is shown 
on a linear contrast scale; the area covered is 200x200 (scale is arbitrary but, for the Solar 
system, it would be in AfJ); and the (single) planet, at r = 30, is shown as the black square 
in Figs. 1 and 2 and as the white square in Fig. 3. 

Fig. 1 illustrates the geometry of dust distribution in different mean motion resonances 
with the planet. One can see that the planet induces dramatic changes in dust distribution, 
compared to the structureless ‘no planet’ case (Sec. 2.1: (i) it creates a highly inhomogeneous 
resonant dust belt that forms near the planet’s orbit, with belt’s width that decreases as j 
increases; (ii) it works to clear from dust the central region within the planetary orbit; and 
(iii) ^ produces an additional highly asymmetric, banana-like empty region or a few such 
regions. 

The above resonant dust belt consists of one or several clumps, which originate as 
follows. In a m : n resonance, each particle experiences n passages through the pericenter 
while the planet makes m revolutions. In the rotating frame, the orbit of the resonant 
particle forms n loops in the vicinity of the planetary orbit. Since in the loop the particle 
moves relative to the planet rather slowly, the particle number density is enhanced there. 

If the amplitude of the azimuthal libration of particle orbits is small compared to the loop 
size, each loop produces one clump. In the opposite case of a large librational amplitude, 
each librating loop creates two clumps located in the two extreme points of libration. 
Thus, in a m : n resonance, the number of clumps equals n at a small and 2 n at a large 
libration. For instance, the 3:1 and 2:1 resonances each produce one clump (Figs. 2b and 
2d), although, owing to libration, it is bifurcated as seen in Figs, la and lb, 2a, 2c, and 2e); 
the 3:2 resonance yields four clumps (Fig. 1c); the 4:3 resonance gives six clumps (Fig. Id); 
the 5:4 resonance produces eight clumps (Fig. le) and the 7:4 resonance - four clumps 
(Fig. If), etc. 

Which of the above resonances is the “strongest’, in terms of the accumulated dust 
mass? The last equation of Sect. 2.2.3 enables us to obtain helpful estimates. For 0 = 0.065, 
one gets Tpn ~ 5 ■ I O' (hereinafter, time is measured in planet’s revolutions, which is 160 
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yrs for a planet at 30 AU). According to our simulations, the lifetime in various resonances 
is as follows: r 4:3 ~ 2 • 10 4 ; ~ 5 • 10 4 ; r 2: t ~ 2.4 • 10 5 ; r 3:1 ~ 8 • 10 s ; r 4:l > 2 • 10 6 , etc. 

If, for simplicity, we put Wj ~ 1 and neglect gravitational scattering of dust grains by the 
planet ( Sj ~ 0), we obtain as reasonable esimates that the dust structures in the resonances 
4:3, 3:2, etc. are comparatively small (Mj/MbZl, whereas resonances 2:1, 3:1, etc. are very 
strong: M } / Mj, 3> 1. Since the probability of capture into the resonance j : 1 decreases for 
larger j , the dust structures associated with the resonances 2:1 and 3:1 seem to happen 
rather often. 


3.2. Dependence of Structure Upon Planetary Mass and Eccentricity 

Given the resonance number, it interesting to examine how the planetary mass 
influences the distribution of dust. For the 2:1 resonance, Figs. 2 a,b and their comparison 
with Fig. lb reveal that as m p i increases, the banana-like empty region increases as well, 
whereas the angle at which the bifurcated substructure is seen from the star decreases: 
virtually one clump might be in the resonance 2:1 with the planet as massive as several mj 

(Fig. 2b). 

While taking into account all the resonances, influence of the planetary mass upon 
the dust structure is more complicated. Qualitatively, at m p //M* — > 0, it is very hard to 
capture dust particles into resonances. At ro p / small but finite, the captured particles are 
in a resonant ring, and if m p i is as small as Earth the planet might be in a resonant lock 
with a trailing clump (Dermott et al. 1994, Reach et al. 1995, TPF 1999). One of the 
important new results of the present study is that the location of the density maximum in 
the wake, i.e. in the trailing zone of the planet, is only valid if the planetary mass is small 
enough; for masses larger, say, than that of Neptune, the density maximum, according 
to our simulations, is located in the leading zone (see Figs. 1 a,b,c,d and Fig. 2). The 
explanation and implications of this phenomenon are given elsewhere (GOM 1999). Here, 
we only notice that (i) depending upon resonance number, the clump(s) may be positioned 
differently relative to the planet - the latter is located either inside the zone of librational 
motion of the loop (Fig. le) or even within the loop - (Fig. If), and (ii) the leading clump 
becomes to be solitary one at the planetary mass as large as a few mj (see Figs. 2b and 2d). 

When the planetary mass is as large as m « 5 mj, it results in dust re-distribution 
in the circumstellar disk shown in Figs. 2 b,c,d. As is seen, the size of the clearing zone is 
appreciably larger, while the zone of libration is much smaller, so that the presence of a 
single clump becomes to be typical. 

Finally, at m pi > 10 mj or so dust again cannot be captured into a resonance such as 
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2:1 due to the overlapping of resonances (e.g. Wisdom 1980). 

Important common features of the models shown in Figs. 1 and 2 include a clearing 
within the orbit of the planet and a banana-like empty region near the very planet (in 
its trailing zone). Dense resonant structures are characterized by richness in features and 
asymmetry. Figs. 1 and 2 give a clear indication of how the dust distribution in the disk 
changes as each of the parameters involved varies, influencing the most essential resonances 
and the number of dust particles captured into them. Asymmetry of the dust cloud 
changes its character for different values of (3 and is larger for smaller (3. Therefore, if 
dust in circumstellar disks is stratified in particle’s size, images of dusty disks should look 
differently, say, in IR and sub-mm wavelengths. 

If the sources of dust are inner, gravitationally scattered comets (like those of 
Jupiter-family in the Solar system), their dust distribution differs drastically from that 
produced by external sources (kuiperoids). Using a stationary distribution of the simulated 
population of Jupiter-family comets computed by us elsewhere (OGT 1998), we can easily 
find the surface density of dust produced by this population . Since the pericenters of 
gravitationally scattered comets are usually located inside the planetary orbit , the dust 
produced by them easily penetrates deeply inside the system. This dust forms a disk-like 
cloud (rather then a belt, as in the case of outer, kuiperoidal sources), with the density 
of dust increasing steeply inward (GOMT 1997). The massive planet makes such a disk 
asymmetric: for instance, Jupiter’s influence results in an arc along the resonant orbit 1:1. 

4. A LINK TO OBSERVATIONS: IMAGERY OF CIRCUMSTELLAR 

DISKS 

There are several points of comparison between our models and observations. The 
most direct comparison is of course the morphology or apparent structure of the disk. 
Besides, our models have implications for interpreting the IR to sub-mm spectral energy 
distributions (SED’s) of circumstellar disks as observed by ISO or SCUBA with the JCM 
telescope, here, we consider some observational signatures of planets embedded in disks. 

Figs. 3a, b show the thermal emission from the simulated disks (seen face-on) with an 
embedded planet, which are to be compared with the available observational data on d Eri 
and Vega shown in Figs. 3c, d (the parent stars are seen almost pole-on). There is a rather 
close resemblance between the observed and simulated disks (the central cavity may be 
masked by emission of dust from non-kuiperodal sources, such as comets). A smaller mass 
induces an asymmetric arc, with two clumps at its edges (Fig. 3a). A similar asymmetric 
.structures has been revealed in sub-mm imagery of the d Eri disk shown in Fig. 3c (Greaves 
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et al. 1998). A more massive planet makes the banana-like region more extended as well 
as merges two clumps into a single one (Fig. 3b). Such a planet may be responsible for 
an asymmetric feature in the simulated disk shown in Fig. 3b, which is reminiscent of the 
observed asymmetric north-east clump near Vega shown in Fig. 3d (Holland et al. 1998). 

Using our modeling, we find that Vega may have a massive planet at a distance of 
50 — 60 AU, and t Eri may have a less massive planet at a similar distance of 55 - 65 
AU. This conclusion is testable: Each resonant feature is stationary in the reference frame 
co-rotating with the planet, but it is not so for the observer at Earth. Therefore, if our 
interpretation of asymmetric clumps as a dynamical resonant structure is correct, the above 
asymmetric feature revolves around the star with an angular velocity of 1.2 — 1.6 deg/year 
(Vega) and 0.6 — 0.8 deg/year (e Eri) - a prediction that can be tested within several 
years. If confirmed, the proposed interpretation of the structure in Vega- and t Eri- like 
circumstellar disks seen face-on would make it possible not just to reveal the embedded 
planet and determine its semimajor axis, but also to constrain its other basic parameters, 
such as mass and eccentricity, and even to pinpoint the position of the planet (see Fig. 3a). 

Although our simulations with just one planet qualitatively reproduce fairly well the 
observed structure of dusty disks in c Eri and Vega, the quantitative details need further 
work. For instance, the present estimates for Vega’s mass, 10 mj, and e Eri’s mass, 0.8 mj, 
(where mj is Jupiter’s mass) should be rather considered as upper limits. Indeed, the 
simulated structure varies with m p ///3 rather than with m p i or (3 alone, and the actual 
value of [3 and therefore m pt might be somewhat exaggerated. A more realistic e also tends 
to diminish ro p /. Collisions of dust particles, which act to populate as well as depopulate 
resonances, is an other factor which needs to be carefully accounted for. We expect that 
further detailed modelling of circumstellar disks in e Eri and Vega would make possible to 
disentangle the basic parameters of the embedded planets. 

In this Letter , we do not consider disks seen edge-on in scattered light, such as the (3 
Pic disc revealed by STIS optical observations (Heap et al. 1999). Neither we discuss here 
dust rings, like HR 4796A or HD 141569, which are addressed elsewhere. 

5. CONCLUSIONS 


As we demonstrated above, a central ‘hole’ void of gas and dust as well as highly 
asymmetric structures in circumstellar disks can serve as indicators of planet(s) embedded 
in the disk. The major goal of the present Letter is to explore how the presence of the 
planet(s) impacts the disk, what particular resonances are responsible for the specific 
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asymmetric features in the dust, and how these features, shape of the disk, etc. would 
enable us to derive the major parameters of the yet invisible planetary system. 

Numerical simulations described above lead to the following major conclusions: 

1. The outermost planet in an exo-planetary system, which mass is > 5 • 10 -5 A/*, can 
produce: (i) central cavity void of dust; (ii) a trailing (sometimes leading) cavity; and (iii) 
an asymmetric resonant dust belt consisting of arc(s) and one, two, or more clumps. 

2. Morphology of these belts and cavities (size, asymmetry, number of clumps, and 
their pattern) depends upon (3 and particular resonances involved and is a strong function 
of the planetary mass, eccentricity of its orbit, and location of dust sources. 

3. Asymmetry of resonant belts increases with mass of the planet. Amplitude of 
libration decreases with the planetary mass and j3. 

4. The crucial test of the above picture would be discovery of revolution of the 
resonant asymmetric structure around the star. For circumstellar disks in e Eri and Vega 
the asymmetric design is expected to revolve by, respectively, (0.6 — 0.8)° and (1.2 — 1.6)° 
annually. 

Numerical simulations described above present a method for determination of planetary 
parameters using a visible morphology of an outher part of circumstellar disk. We have 
demonstrated above that a very rich structure of the dusty disk with an embedded planetary 
system and its dependence upon the parameters of the planet and dust offer a possibility 
to reveal the presence of the planet(s) and to disentangle the planetary orbital parameters 
from dust parameters. Moreover, our high-resolution simulations ensure an opportunity to 
evaluate the major planetary characteristics, such as mass, semimajor axis and eccentricity. 

The method of determining the planetary parameters using the shape and structure 
of the circumstellar disk is complementary to the current approaches of finding exo-solar 
planets using precise Doppler measurements, microlensing, and transits. Those approaches 
operate for Jupiter-like masses located close to the parent star so that the orbital period 
usually does not exceed several years. Our method yields information on planets including 
the outermost ones located at tens or hundreds of AU (i.e. with periods of several hundred 
years) and having variety of masses, both smaller and larger than Jupiter. 

In sum, our modeling offers the ability to determine the major orbital parameters and 
masses of planets in dusty disks. As these disks are common, many planetary systems 
may be found. Moreover, since the dust disks are often so bright that they prevent direct 
detection of planets with imaging at visible, [R, or sub-mm wavelengths, our technique 
may actually be helpful in finding more planets than other methods do. ll is well worth 
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the efforts to develop the modeling and interpretation for observations already available or 
obtainable in the near future. 
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TABLE 1 


Comparison of Our Approach With That of RSSS 1 



imber of particles in the grid 


RSSS 

8192 


Our approach 
l 0 i°-u 


-ty of resonances accounted for 

Relevance to planetary disks 

) Computationals errors 
md integrator stability 

Constraints on time scale 
of evolution, r evo / 

6 ) Computing facility, 

CPU time 

No. of different variants 
ange of planet’s mass 

iparability with observations 


a few low-order resonances 

particles of the same age 

4th order Runge-Kutta 
integrator with single precision 

T evoi < 10 4 revolutions of the planet 4 

Connection machine CM 2/8k, 

6 hrs 

9 

10 -3 0.15 m j 

no disk emission considered 


all mean motion resonances 

particles of all ages 

implicit 2nd order Potter’s 
method with double precision 1 2 3 ' : 

Tevoi < 2 • 10 6 revolutions of the pla 

several PCs 
a few to a few tens hrs 

40 

0.05 m j^rript t 10 m j 
disk emission considered 


1 Roques, Scholl, Si.cardy k Smith (1994) 

2 Taidakova (1997) 

3 Taidakova k Gorkavyi (1997) 

RSSS neglected 20% of long-living resonant particles 
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